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In effective models of loop quantum cosmology, the holonomy corrections are associated with 
deformations of space-time symmetries. The most evident manifestation of the deformations is the 
emergence of an Euclidean phase accompanying the non-singular bouncing dynamics of the scale 
factor. In this article, we compute the power spectrum of scalar perturbations generated in this 
model, with a massive scalar field as the matter content. Instantaneous and adiabatic vacuum- 
type initial conditions for scalar perturbations are imposed in the contracting phase. The evolution 
through the Euclidean region is calculated based on the extrapolation of the time direction pointed 
by the vectors normal to the Cauchy hypersurface in the Lorentzian domains. The obtained power 
spectrum is characterized by a suppression in the IR regime and oscillations in the intermediate 
energy range. Furthermore, the speculative extension of the analysis in the UV reveals a specific 
rise of the power leading to results incompatible with data. 

PACS numbers: 98.80.Qc, 98.80.Jk 


I. INTRODUCTION 

Loop quantum gravity (LQG) is a simple, consis¬ 
tent, non-perturbative and background-independent 
quantization of general relativity. It uses Ashtekar 
variables, namely the SU(2)-valued connections and 
the conjugate densitized triads. The quantization is 
obtained through holonomies of the connections and 
fluxes of the densitized triads. No heavy hypothesis is 
required. Introductions can be found in Refs. PQ. 

Loop quantum cosmology (LQC) is an application of 
LQG-inspired quantization methods to a gravitational 
system with cosmological symmetries. In LQC, the 
big bang is generically replaced by a big bounce due 
to repulsive quantum geometrical effects when the 
density approaches the Planck density and interesting 
predictions can be made about the duration of inflation 
when a given matter content is assumed. It is, however, 
important to underline that LQC has not yet been 
rigorously derived from LQG and remains an attempt 
to use LQG-like methods in the cosmological sector. 
Introductions can be found in Refs. [a ej. 
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The confrontation of LQG with available empirical 
data is crucial in order to check the physical validity of 
this approach to quantum gravity. The most promis¬ 
ing option in this direction is currently given by the 
exploration of the cosmological sector of LQG. The 
present state of advancement, however, does not allow 
for a derivation of the cosmological dynamics directly 
from the full theory. Because of this, LQC models are 
considered to fill the existing gap. These models suffer 
from quantum ambiguities, which are believed to be 
fixed by the cosmological dynamics regained from LQG. 

This study is based on the effective LQC dynamics, 
which allow to address various cosmological issues. In 
particular, numerous studies have been devoted to the 
computation of tensor power spectra and their signifi¬ 
cance in the light of the future observations (see, e.g., 
Refs. [3H7] ). In this work, we will focus on scalar modes, 
which are more relevant from the observational point 
of view but which are more demanding to deal with at 
the theoretical level because of subtle gauge-invariance 
issues and hypersurface deformation algebra closure 
conditions. 

Two main types of quantum corrections are expected at 
the effective level of LQC. The first one comes from the 
fact that loop quantization is based on holonomies, i.e. 
using exponentials of the connection rather than direct 
connection components. The second type of corrections 
arises for inverse powers of the densitized triad, which, 
when quantized, becomes an operator without zero 



eigenvalue in its discrete spectrum, thus avoiding the 
divergence. As the status of “inverse volume” corrections 
is not fully clear, due to the fiducial volume dependence, 
this work focuses on the holonomy term alone which has 
a major influence on the background equations and is 
better controlled [8]. In this framework, we will consider 
the Euclidean phase predicted by LQC B mg, and 
put, as advocated in Ref. m, initial conditions in the 
remote past of the contracting branch of the universe 
(this choice can be questioned and other proposals have 
been considered mm- 

It is worth noticing, that an alternative attempt re¬ 
garding the cosmological perturbations in LQC have 
recently been presented (see Refs. jH]). In this ap¬ 
proach, quantum fields are considered on a homogeneous 
quantum background, based on the methods developed 
in Ref. 15]. Because the gauge-invariant variables for 
perturbations are fixed to be the classical ones, the 
Euclidean phase characterized by the elliptic nature of 
the equations of motion does not occur. However, the 
consistency of the effective dynamics emerging from this 
formulation remains an open issue. 

The key ingredient of this work is the existence of 
an Euclidean phase around the bounce. A completely 
rigorous study of the physical consequences would 
require a full understanding of quantum field theory 
(QFT) in curved Euclidean spaces and at the junction 
hypersurface with the Lorentzian manifold. This is 
obviously far beyond the scope of this article. Especially 
when taking into account that even on a well behaved 
Lorentzian dynamical space, QFT is not without am¬ 
biguities and many issues remain open. The central 
methodology of this study is to define physical quantities 
in the initial classical Lorentzian space (the contracting 
branch before the bounce), to evolve them in a mathe¬ 
matically rigorous way through the Euclidean zone, and 
to calculate observables in the second Lorentzian phase 
(the expanding branch we live in) where the physical 
meaning is clear again. Although, by definition, there 
is no time anymore in the Euclidean phase, one can 
still rely on general covariance and the corresponding 
structure has a well-defined canonical formulation using 
hypersurface deformations. We don’t claim that this is 
the only way to address this situation. There are clearly 
other possible ways to deal with this speculative new 
phenomenon. However, it seems to be a quite natural 
and reasonable first assumption. In addition, this is the 
methodology that has been used up to now to evaluate 
the tensor spectrum in the deformed algebra approach 
to LQC. It is therefore important to also derive the 
scalar spectrum following the very same methodology, 
at least for a meaningful comparison. 

In the following section, we first remind the basis of the 
deformed algebra approach used in this study. In section 
III, we summarize some important features of the back¬ 
ground dynamics in LQC. The equation of motion for 


scalar perturbations is derived in section IV. In section 
V, different ways of choosing initial conditions for per¬ 
turbations are presented. Section VI is devoted to the 
analysis of the scalar power spectrum. Concluding re¬ 
marks are given in section VII. 

II. DEFORMED ALGEBRA 

In the canonical formulation of general relativity, the 
Hamiltonian is a sum of three constraints, 

H g [N\ N a , N] = - 1 - [ d 3 x (VC, + N a C a + NC) « 0, 
2k Js 

where k = 8ttG , ( N l ,N a ,N ) are Lagrange multipliers, 
Ci is the Gauss constraint, C a is the diffeomorphism con¬ 
straint and C is the scalar constraint. The equality de¬ 
noted as is to be understood as an equality on the 
surface of constraints [i.e. a weak equality). It is conve¬ 
nient to define the corresponding smeared constraints, 

C 1 = G[N i } = ^ J^d 3 x WCi, ( 1 ) 

C 2 = D[N a ] = ±-J^d 3 x N a C a , (2) 

C 3 = S[N] = ljd 3 x NC, (3) 

such that Hg[N\ N a , N] = G[W] + D[N a } + S[N]. 
The Hamiltonian is a total constraint which vanishes 
for all multiplier functions (N l ,N a ,N). The time 
derivative of the Hamiltonian constraint vanishes also 
weakly and therefore the Hamilton equation, / = 

{/, H G [M i , M a , M] }, leads to 

{H G [N i ,N a ,N],H G [M i ,M a ,M]} « 0. (4) 

As the Poisson brackets are linear, the condition Q is 
satisfied if the smeared constraints belong to a first class 
algebra, 

{Ci,Cj} = f K ij(Ai,Ef)C Ki (5) 

where the f K IJ (Ai,E^) are structure functions which 

depend on the Ashtekar variables (A J b ,Ef). The alge¬ 
bra closure is fulfilled at the classical level due to general 
covariance. The algebra must also be closed at the quan¬ 
tum level. Otherwise the system might escape from the 
surface of constraints, leading to an unphysical behavior. 
In addition, as shown in Ref. I6j , the algebra of effective 
quantum constraints should be strongly closed (that is, 
off shell closure must be considered). This means that the 
relation ([5| should hold in the whole kinematical phase 
space, and not only on the surface of constraints (cor¬ 
responding to on shell closure). When the constraints 
are quantum-modified by the holonomy corrections, the 
resulting Poisson algebra might not be closed, 

= f K u(Ai, E?)C% + Au , (6) 
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where Aij stands for the anomaly term which can ap¬ 
pear due to the quantum modifications and the super¬ 
script ‘Q’ indicates that the constraints are quantum cor¬ 
rected. The consistency (closure of the algebra) requires 
that all Aij should vanish. Remarkably, the conditions, 
Aij = 0, lead to restrictions on the form of the quantum 
corrections and determine them uniquely under natural 
assumptions. 

The issue of anomaly freedom for the algebra of cosmo¬ 
logical perturbations was extensively studied for inverse- 
triad corrections. It was demonstrated that this require¬ 
ment can be fulfilled at first order in perturbations for 
scalar mm , vector m and tensor perturbations m- 
Predictions for the power spectrum of cosmological per¬ 
turbations were performed [21] , leading to constraints on 
some parameters of the model by the use of observations 
of the cosmic microwave background radiation (CMB) 
[22] . It was also considered for holonomy corrections and 
vector modes in Ref. [23] and for scalar modes in Ref. m 
(the full analysis with inverse-triad and holonomy terms 
was performed in Ref. [24] ). It was shown in Ref. ,25] that 
there exists a single modification of the algebra structure 
that works for all kinds of modes, thus emphasizing the 
consistency of the theory. It is also important to un¬ 
derline that the matter content plays a role in removing 
degeneracies. Even if the calculations carried out in the 
above-mentioned articles are quite laborious, the guiding 
idea behind is very simple. Each time a k factor, de¬ 
fined as the mean value of the Ashtekar connection A l a , 
appears, it is replaced by 


k —> 


am(nfjrfk) 

nfi'Y 


(7) 


where n is some unknown integer and p is the coordinate 
size of a loop. The full perturbations have to be cal¬ 
culated up to the desired order, the Poisson brackets are 
then explicitly calculated and the anomalies are cancelled 
by counter-terms required to vanish in the classical limit. 
The neat result is that the algebra of effective constraints 
is deformed with respect to its classical counterpart. It 
takes the following form: 

{D[M a ], D[N a ]} = D[M b d b N a - N b d b M a ], 

{D[M a ],S Q [N]} = S Q [M a d b N - Nd a M a ], 

{S q [M], S q [N]} = SID [q ab {Md b N - Nd b M)\ , 

where Q is the deformation factor that plays a crucial 
role in the following. It is given by f2 = 1 — 2 p/p c where 
p is the density of the Universe and p c is the critical 
density expected to be close to the Planck density. In 
Lorentzian General Relativity fl = 1. When Q < 0 
the structure of space-time becomes Euclidean. (Strictly 
speaking space-time is Lorentzian or Euclidean only if 
O = ±1 but the most important properties regarding 
physical consequences, namely the existence of a causal 
structure and the general behavior of the solutions for 
wave equations, only depend on the sign of Q and not on 
its precise value [13] ■ It therefore makes sense to speak 


of Lorentzian or Euclidean phases.) 

Interestingly, this conclusion has strong links with results 
often postulated (for technical reasons, notably a better 
behavior of path integrals) in quantum cosmology, but it 
appears here as a real dynamical prediction of the theory. 
In standard quantum cosmology one usually deals with 
an amplitude 

<02, *21 <£i ,ti>= y^[0]e /[01 , (8) 

where 7[0] is the action of the field configuration 
and d[0] is a measure on the space of field configura¬ 
tions. The integrand of Eq.(8]) has a rapidly oscillating 
phase, and the path integral, in general, does not con¬ 
verge. This is why the time is rotated clockwise by 7t/2 
so that 7[0] —> 7[0] = — il[<f>}. The integrand in the 
resulting Euclidean path integral is now exponentially 
damped, and the integral generically converges. Then, 
one can analytically continue the amplitude in the com¬ 
plex f-plane back to real values. Importantly, a quan¬ 
tum field theory machinery has been developed in this 
framework and the interested reader can, e.g., consider 
Ref- [2?]■ Of course, there are also links with the Hartle- 
Hawking proposal [25]. But the Euclidean phase appears 
in the model considered in this article in a fundamentally 
dynamical way since the Poisson bracket between Hamil¬ 
tonian constraints varies continuously from a positive to a 
negative expression. The “spirit” of the Hartle-Hawking 
proposal, translated in the framework considered here, 
has been studied in m ■ Importantly, the appearance 
of an Euclidean phase was also independently derived 
from another approach to LQC in Ref. [2Sj- Many quan¬ 
tum gravity approaches seem to predict the existence of 
a silent surface (fi = 0) where light cones are completely 
squeezed, on each “side” of the Euclidean phase. This is 
also a clear realization of the BKL conjecture (see, e.g. 
Ref. [29]). Arguments are given in Ref. [8j showing that 
the change from a hyperbolic to an elliptic type of equa¬ 
tions in LQC should be understood as a true change of 
signature (that has been missed before because homoge¬ 
neous models cannot probe it) and not just a tachyonic 
instability. It should also be emphasized that the de¬ 
formed algebra approach is grounded in avoiding gauge 
issues. Many approaches make a gauge fixing. In most 
cases, gauge fixing before quantization is known to be 
harmless, but the situation is different in general relativ¬ 
ity. The constraints we are considering are much more 
complicated functions than, for example, the Gauss con¬ 
straint of Yang- Mills theories: it is therefore likely that 
the constraints receive significant quantum corrections. 
If the constraints are quantum corrected, the gauge trans¬ 
formations they generate are, as we have shown, not 
of the classical form. Gauge fixing before quantization 
might then be inconsistent because one would fix the 
gauge according to transformations which subsequently 
will be modified. In addition, in the present case, the 
dynamics is part of the gauge system. A consistent the¬ 
ory must therefore quantize gauge transformations and 
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the dynamics at once. It is not correct to fix one part 
(the gauge) in order to derive the second part (the dy¬ 
namics) in an unrestricted way. The subtle consistency 
conditions associated with the covariance of general rel¬ 
ativity are encoded in the first class nature of its system 
of constraints. Here, great care is taken in not breaking 
this consistency. 

The equations of motion derived in this framework are 
still covariant under the deformed algebra replacing clas¬ 
sical coordinate transformations. The corresponding 
quantum space-time structure is obviously not Rieman- 
nian (there is no line element in the usual sense), but has 
a well-defined canonical formulation using hypersurface 
deformations. 


III. BACKGROUND EVOLUTION 


At the background level, the change of signature can¬ 
not be probed/detected. This is obvious for two rea¬ 
sons. First, because the relative sign between tempo¬ 
ral and spatial derivatives cannot be identified. Second, 
because the Poisson bracket between Hamiltonian con¬ 
straints then trivially vanishes. The evolution of the cos¬ 
mological background is studied at the effective level with 
holonomy corrections. The background geometry is de¬ 
scribed by the homogeneous, isotropic and flat configu¬ 
ration parametrized by the scale factor a. The dynamics 
of the background is governed by the quantum-corrected 
Friedmann equation 



(9) 


derived in Ref. [3D], where H = a/a is the Hubble rate 
in cosmic time, p is the energy density of the content of 
the universe and p c denotes its maximal value attained 
at the bounce. The dot denotes a derivative w.r.t. cos¬ 
mic time. Obviously, the physical interpretation of the 
background equations makes sense only in the Lorentzian 
phase, that is when p < p c / 2, but technically one can 
still determine the manifold structure for p c / 2 < p < p c , 
that is in the Euclidean phase. An alternative view of 
our approach which solves all interpretation difficulties 
and which leads to the very same result, is to consider 
that there is no real change of signature: the background 
evolves in a standard way (although, of course, according 
to the modified Friedmann equation) and the change of 
sign of the fl factor entering the propagation equation 
of perturbations is just a tachyonic instability, as well 
known to appear, for example, in Gauss-Bonnet gravity 
when the curvature invariant is non-minimally coupled 
with a scalar field. We prefer not to favor this view be¬ 
cause, as detailed in [S], there are hints that the phe¬ 
nomenon is deeper but one is free to see things in this 
way. Planck units are used throughout this article with 
rripi = l/y/G « 1.22-10 19 GeV. Furthermore, we consider 
a single massive scalar field, 0, with a quadratic poten¬ 
tial V = m 2 cj) 2 / 2, as the matter content of the Universe. 


This choice is made for simplicity. It allows easy compar¬ 
isons with other works and generates a phase of slow-roll 
inflation. Even if this potential is not favored by cur¬ 
rent observational data [31], it still serves as a valuable 
toy model for studying the phase of inflation in different 
frameworks. Moreover, taking into account more subtle 
effects, e.g. the quantum gravitational corrections con¬ 
sidered here, might improve the status of the quadratic 
potential in the light of the observational data. 

Splitting the field (j> = <j> + 8<t> into a background part, </>, 
and a perturbed part, 5<p, the Klein-Gordon equation for 
the background reads 

4> + 3H$ + m 2 ^ = 0. (10) 

A first analysis of this model has already been studied in 
Ref. [32]; a detailed analysis of the background equations 
can be found in Ref. [33]. Here, we only summarize the 
main features of the background dynamics. The field evo¬ 
lution can be characterized by two dynamical parameters 
[6], the potential energy parameter, x, and the kinetic en¬ 
ergy parameter, y, defined as 


rruj) 4> 

:= 72 pT’ y:= vW: 


( 11 ) 


Then the total energy density can be written as p = 
p c {x 2 + y 2 ). Eqs. § and (jTo]) can then be recast as 

( H = —np c y 2 (l — 2a: 2 — 2y 2 ) , 

< x = my , ( 12 ) 

[y =-3Hy-mx, 

showing that there are two timescales involved in this 
system: one is given by 1 /m and corresponds to the clas¬ 
sical evolution of the field, the other one is 1 /y/3np c and 
corresponds to the quantum regime of the evolution. The 
ratio of these two timescales is 


r := 


m 

V3 Kp c 


(13) 


According to standard assumptions of slow-roll infla¬ 
tion with a quadratic potential, the value of the mass 
m~ 1.2 x 10~ 6 m P1 is preferred in the light of the obser¬ 
vational data from the Planck satellite (see Ref. j3T]). 
The critical energy density at the bounce is given by 
p c = 0.41 TOpj, which is exactly the upper bound of the 
spectrum of the energy density operator [3]. These val¬ 
ues lead to r ~ 2 • 10~ 7 . Hence, one can safely assume 
that r<C 1, ensuring that the evolution splits into three 
phases: (i) a classical pre-bounce contracting phase, (ii) 
the bouncing phase and (iii) a classical expanding phase 
after the bounce (slow-roll inflation), see Ref. [33] for de¬ 
tails. Initial conditions, {a 0 , x 0 , y 0 }, are set in the remote 
past of the contracting phase when the energy density is 
very small compared to the critical energy density, i.e. 



(14) 
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It is convenient to use polar coordinates for the potential 
and kinetic energy parameters 


x(t) = ./—sin (mt + 9 0 ), 

(15) 

V dc 

y(t) = \ — cos (mt + 6 0 ), 

(16) 

V Po 


where 0 o is the initial phase between the share of poten¬ 
tial energy and kinetic energy. In order to select different 
background evolutions independently of the small oscilla¬ 
tory behavior of the solutions, the following parametriza- 
tion shall be used: 


/Po = T / _ sin (20 o ) 


a 


4a 


(17) 


where a is a number large enough such that (14) holds. 


To each phase, 9 a , corresponds a specific value of the po¬ 
tential energy parameter at the bounce x B . As shown in 
Ref. [M], for a mass m = 1.21 x 10 _6 to p1 , the favored 
value for x B is 3.55 x 10~ 6 . This solution for the back¬ 
ground dynamics features only a tiny amount of deflation 
before the bounce as shown in Fig. [lj In general, we will 
chose the normalization of the scale factor at the bounce 
as a B = 1. The plots and spectra are presented as func¬ 
tions of the number of e-folds N := ± ln(a/a B ), that have 
to elapse until the bounce (negatively valued) and that 
have elapsed after the bounce (positively valued) respec¬ 
tively. 





FIG. 1. Evolution of the scalar field as a function of the 
number of e-folds N := ±ln a/a B , with m = 1.2 x 10~ 6 m P i. 
The zero on the horizontal axis corresponds to the bounce 
when a B = 1. This solution is such that x B = 3.55 x 10 -6 
(obtained with a = 177r/4 + 1 and 9 0 = 5.11). The evolution 
is stopped at the end of inflation when cj) = l/\/47r m P i. 


IV. EQUATION OF MOTION FOR SCALAR 
MODES 


form of the Hamiltonian constraint. In Ref. [2S], the 
gravitational part of the Hamiltonian constraint has been 
analyzed and reads (up to quadratic order) 


H[N\ = 



N(H (0) +H {2) ) + SNH ( ' 1) 


(18) 


where 

2 kH W 

2 kM (1) 


2k n {2) 


-&Vpk 2 , (19) 

-4 y/§SKi - ^SEj + A. d j d c 6E °, (20) 

-2^=6 K'JE? 

Vp 


+ Vp(S b i 6K' a 6»6Ki - StSK^SKi) 

+ —(SISE^ISE^ 28i6E?Si6E$) 

2 

+ in'f da (SEfdiSE?) 

P 2 

+ {d c 6E?)(d d 6E$). (21) 

P 2 


Here, p is the mean value of the densitized triad Ef and 
k was defined in Eq. ([7]). The term Z'^f 3 depends on the 
kind of modes considered (scalar, vector or tensor): 


{ SabS^S 04 for tensor modes, 

0 for vector modes, (22) 

— for vector modes. 


Lastly, Y^ 1 is a complicated expression whose form is 
not relevant here. Based on this, the holonomy quan¬ 
tum corrections can be accounted for and the Mukhanov- 
Sasaki equation of motion for gauge-invariant perturba¬ 
tions can be calculated m- In conformal time it is given 
by 

z " 

«"-SlV 2 i) s -A» s = 0, (23) 

z s 

with 

V S ■■= Vp (&<!>+ and := Vpjf- (24) 


The variable $ denotes the gauge invariant Bardeen 
potential taking into account the metric perturbations, 
whereas tf> represents the massive scalar field. is the 
conformal Hubble parameter. The Mukhanov equation 


of motion (231 reduces to the classical equation when 


fl —> 1. Note that for FLRW cosmologies, in conformal 
time, yjp = a. On the quantum-modified background 
discussed in the previous section, we can evaluate the 
evolution of the Mukhanov variable v s . For simplicity 
we will omit the index ‘S’ in the following, assuming that 
it is clear that v denotes the scalar perturbation variable. 
Using the Fourier space decomposition of the u(x, rj) field, 


;(x,r?) = 


d 3 k 

(2tt) 3 /2 




The equation of motion for scalar modes in the de¬ 
formed algebra approach is derived from the particular 
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(25) 









one gets a set of ordinary differential equations for the 
Fourier components Uk- Due to the isotropy of space, the 
k-vector of i’k might be simplified to the absolute value 
dependence v k , where k := \/k ■ k. The v k function is 
called a mode function. Instead of using the conformal 
time dependence, it is often (due to technical reasons) 
convenient to switch to cosmic time t in the numerical 
computations. With t = J a- dr), the Mukhanov equation 
of motion reads 


v k + Hv k + fk\t)v k = 0, 


(26) 


with z = a and 


/T(t) 


(27) 


being the effective frequency term. In order to derive 
the primordial power spectrum after inflation one would 


have to solve Eq. (26) for every mode k for all times from 
f init until f end where t init is the initial starting point, set in 
the remote past as we will see later, and t end denotes the 
time at the end of the inflationary phase. This requires a 
numerical integration. However, the Mukhanov variable 
v. cannot be used for the whole integration because of a 
non-physical singularity occurring at the bounce. Let us 
describe how to bypass this difficulty by using the change 
of variable. We introduce hk := Vk/a for every k , so that 
(26) becomes 


hk + 3 Hhk + fl ( t)hk — 0 


(28) 


with 


fk\t) :=n^+m 2 + m 2 Kfi^-2 ( ^ j +^. (29) 


In the numerical computation this second order differ¬ 
ential equation is replaced by the following first order 
system: 


h k = (1 /a)g k , 

g k = -2 Hg k + f k \t)a h k - 


(30) 


The numerical integration of (30) is performed for t £ 
where is before the bounce. Since the 


effective frequency terms, (29) or (27), depend on inverse 


powers of H , the differential equations have a generic sin¬ 
gularity at the bounce, when the Hubble parameter van¬ 
ishes. Nonetheless, this singularity is not physical, which 
can be seen by analyzing the physical scalar curvature, 


v 

iz := 

z 


(31) 


Using this variable, one can rewrite Mukhanov’s equation 
of motion in Fourier space as 


Hk ~ 1 3il + 2 m 2 t + 2^ ) h k + n^TZk = 0. (32) 

V rhU a z 


When approaching the bounce, H tends to zero and fl 
to minus one. It should be noticed that this equation is 
mathematically well behaved, even when Q < 0, that 
is in the Euclidean-like phase. Obviously, in the Eu¬ 
clidean phase, k loses its usual interpretation and the 
Fourier transformation has no intuitive physical sense. 
This is, however, already true on a Lorentzian curved 
manifold: the space/ time splitting - or positive/ nega¬ 
tive frequency splitting - is, in general, ill defined: one 
usually relies on static boundaries that allow a physi¬ 
cal choice. The “bet” of this study is fundamentally 
the same: quantities are defined in the remote past of 
the contracting branch when the signature is the usual 
one, and when quantum effects are negligible. Then, 
they are mathematically propagated through the bounce 
where they lose a clear meaning. Finally, the observables 
are computed in the expanding Lorentzian phase when 
physics is again under control. This is a questionable 
approach, but, in our opinion, a reasonable one at this 
stage. 

During the bouncing phase, </> ~ a/ 2p c mcj), so that 

the equation of motion reduces to 

H ■ k 2 

h k — 2—hk - ~TZk — 0. (33) 

H a z 

According to the analytical expression for (j) and <f> around 
the bounce developed in Ref. [33], this gives 

h k - hk - k 2 TZk — 0, (34) 

t t B 


when we restrict ourselves to the first order in (t — t B ). 
The space of solutions to this differential equation is 
spanned by the two independent functions, 

1Z / x) = [sinh (k(t — t B )) — k(t — t B ) cosh ( k(t — t B ))], 

7= [cosh ( k(t — t B )) — k(t — t B ) sinh ( k(t — t B ))] • 


These solutions show an obviously regular behavior at 
the bounce. But it should be noticed that (32) runs into 


trouble away from the bounce due to the time deriva¬ 
tive of the potential (f) appearing in the denominator of 
the friction term. During the classical contracting and 
expanding phases, </ oscillates around a null value, caus¬ 
ing the break-down of the numerical integration of the 
differential equation (32) of 7 Z. For this reason one has 


to switch twice between Eq. (30) and Eq. (32) during 
the numerical computations. For t £ [t iIlit ,t h ^, n ] and 
t £ [tn^hhen d]> where t h ^ u ^t B < t^ h , the differen¬ 


tial equation for h, namely (30), must be used. Whereas 
for t £ [t h -> n ,t u ^ h ], it is the equation for 7 Z, Eq. (32), 


that has to be integrated. The exact choice of the transi¬ 
tion points is irrelevant as long as they do not approach 
one of the singularity points. 
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V. INITIAL CONDITIONS 

The considered equations of motion for the scalar per¬ 
turbations (Eq. |23j ) or Eq. (26)) might be considered 
- at the effective level as the quantum ones. This is 
due to the presence of the factor S7, being a result of 
the quantum gravitational effects. The quantum effects 
taken into account here are, however, only those which 
modify the background degrees of freedom. The inhomo¬ 
geneous degrees might be (and are), treated classically 
in the perturbative regime under consideration. In or¬ 
der to see it, let us consider the extrinsic curvature K l a 
which is exponentiated to the form of a holonomy opera¬ 
tor in the quantum theory. In the perturbative treatment 
we have K l a = k8\ + SK *, together with the condition 
\8K l a \/k -C 1. Path integration of K l a leads to a factor of 
the form 7 jlk for the homogeneous contribution, which 
is of the order of unity in vicinity of the bounce. Full ex¬ 
ponentiation of the background contribution to K' a must 
be, therefore, kept over the evolution through the bounce. 
However, this is not necessary for sufficiently small per¬ 
turbations, for which the condition 7 /l|( 5 /\*| <C 1 might 
be satisfied even at the bounce. This allows for the ex¬ 
pansion of the holonomy up to a linear contribution in 
SK l a and treating the perturbative degrees of freedom in 
a classical manner. 

As the phase space of the perturbative degrees of free¬ 
dom is approximated by the classical one, the canonical 
quantization procedure for the modes might be applied. 
The canonical quantization is an approximation, which 
is valid only for sufficiently small amplitudes of 8K l a and 
sufficiently large wavelengths (roughly greater than the 
Planck length) in the mode expansion. If the conditions 
are satisfied, the Fourier mode v k (rj) can be promoted to 
be an operator, such that in the Heisenberg picture 

Dc (v) = Vk(v)a* + vUv)aL k , (35) 

where 17 ( 7 ) are the so-called mode functions satisfying 
the classical equation (26). The and a k are the cre¬ 


ation and annihilation operators respectively, satisfying 
[a k ,dq] = <^ 3) (k — q)- Using this commutation relation, 
one may show that the following condition 


dv 


k * dv k 


Vk Af~ V ^ =tt 


(36) 


called Wronskian condition, has to be satisfied in order 
to preserve the standard canonical structure. 

Based on the above, the two-point correlation function 
for the scalar curvature field 7?.(x, rj) in the vacuum state 
| 0 ) is given by 


(0|7t(x, T])il(y, r))\0) = 


r °° dk T> a \ sinkr 

-r V s(k,r ])——, (37) 

l AC A zr 


where r = |x — y| and the scalar power spectrum reads 

k 3 \vk{v) I 2 


The power spectrum carries all statistical information 
about the Gaussian scalar curvature field under consid¬ 
eration (non-linear effects are neglected in our analysis) 
and its determination for the LQC model discussed in 
the previous sections is a main goal of this study. 

The equations of motion for the scalar mode functions 
17 ( 77 ) can be solved numerically, following the procedure 
presented in the previous section. For this purpose ini¬ 
tial conditions for perturbations have to be set for every 
wavenumber k. In standard cosmology, it is common to 
set Cauchy initial conditions at some moment in time 
after the big bang singularity. In the present work, we 
set initial conditions in the pre-bounce phase. This is 
the natural choice if the bounce is a phenomenon to be 
really understood as resulting from a causal evolution of 
the Universe, with time flowing in a unique direction. In 
addition, in the remote past of the contracting branch, 
the Universe is classical and quantum effects do not play 
any important role. This seems both technically more 
convenient (since the quantum dominated region still rep¬ 
resents a quite unknown held of physics) and physically 
better motivated. In particular, as discussed in the pre¬ 
vious sections, it has been shown in Ref. m that the 
geometry of the universe in its quantum stage (p > p c / 2 ) 
might become Euclidean instead of being Lorentzian (the 
very notion of time obviously looses here its meaning). 
The physical consequences are still not perfectly well un¬ 
derstood and setting initial conditions in the Euclidean 
phase would be the worst possible choice: we therefore 
focus on the classical contracting phase. Note that an¬ 
other proposal, studied in Refs. tuna, is to set initial 
conditions at the surface of silence (or in a “hybrid” way 
as advocated by a careful study of the Tricomi problem). 
We do not consider this hypothesis here. 

The most simple and natural way to set initial condi¬ 
tions for a quantum oscillator is provided by the vac¬ 
uum state | 0 ) fc for every mode k at some given moment 
in time. This sets a clear Cauchy initial value prob¬ 
lem. However, it is well known that the notion of vac¬ 
uum in an arbitrary curved spacetime is ambiguous since 
the definition of the usual “instantaneous vacuum” is 
based on plane waves satisfying the differential equation 
of an harmonic oscillator with wavenumber k. In the 
case of scalar perturbations, as the effective frequency 
term depends non-trivially on time, it is more appro¬ 
priate to use the Wentzel-Kramers-Brillouin (WKB) ap¬ 
proximation and the so-called “adiabatic vacuum”. In 
contrast with the ordinary instantaneous vacuum state, 
well known from quantum held theory, the adiabatic vac¬ 
uum does not have to satisfy the differential equation of 
an harmonic oscillator in some limit. The instantaneous- 
vacuum state is recovered as the first term of the WKB 
expansion. 

For scalar perturbations in LQC the equation of motion 
is given by Eq. (231. In conformal time, this equation 


'Ps(k, rj) := 


2-7T 2 


(38) 


resembles the differential equation of an oscillator, 


v'k + k l „(f?H = 0, 


(39) 
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with a time-dependent wave number 


fceffO?) : = \!^{vW - —(??). 


(40) 


Recall that we set initial conditions in the remote past 
where 17 ~ 1 is almost constant. Thus the 17-factor ac¬ 
tually plays no role when addressing the issue of initial 
conditions. The main idea of the WKB approximation 
is to use the following generic ansatz for the solutions to 
Eq. _ 


Vfc(jy) = ci • e i ^ liT) - WkM + c 2 • e~ i(KttT) Wk{r, \ (41) 

where the values of constants Ci and c 2 are constrained 
according to the Wronskian condition (361. In the 


WKB approximation the functions Wk{jf) are expanded 
in terms of some small parameter (fc eft T)~ 1 where T is the 
minimal time interval for which k e!! , and its time deriva¬ 
tives, start to change substantially ( T l/k e!t ). Then 
the WKB expansion reads 


w k (v) = 


n —0 


Ki C t 


W k , n (r]). 


(42) 


Introducing this ansatz into Eq. (39), one gets the ex¬ 
plicit expressions for the different orders, n, of W k , n - The 
WKB approximation consists in truncating the series af¬ 
ter the first order, leading to the approximated solution 
for the mode functions 


v k (r]) = Cl P i P kett(v)dfj 


VKhJ rj) 


+ 


° 2 r —i f v Kft(fi)dfi 


VhJrjj 


(43) 


Using the Wronskian condition (36), we find |c 2 |~ — 


| ci | = 1/2 as a condition that the free parameters Ci and 
c 2 have to fulfill. The most convenient choice is ci = 0 
and c 2 = l/y/2 which corresponds to a wave propagat¬ 
ing in positive time direction. For this choice the mode 
function reads 


v k (v) = 


1 


\/2 kMrjj 


,—i J v keff(ij)dfi 


(44) 


This represents a suitable choice because the mode func¬ 
tion reduces to the Bunch-Davies vacuum in the UV limit 
(k —> oo). 

In order to check that this approach is valid, we plug this 
solution into Eq. (39) and find that it is actually an exact 
solution to 


- 

"'eff 


3 (fc( ff ) 2 

4 kl t 


1 k" 

2 K 


v k = 0. 


Therefore, the solution (43) is valid as long as 


1C 

2 *£ t 


3 (fc(tf) 2 

4 


« 1, 


(45) 


(46) 


when the effective wavenumber, fc eft , varies slowly. The 
appropriate initial conditions are then 

1 


«fc(?7i„it) = 


V^AW^init) 


dv k 

dr] 


V=ViTiit 


= ~ iK«+- 


1 K 

2 k eff 




(47) 


(48) 


V=Vinit 


where the exponential term can be neglected because it 
contributes only with an arbitrary phase. The initial 
moment has to be chosen such that the WKB condi¬ 
tions are satisfied at this particular moment of time for 
all modes k. By analyzing fc eff (r/) and its time derivatives, 
we can find an appropriate ry init in the remote past. For 
the numerical computations, the choice of rj init is there¬ 
fore arbitrary as long as the condition (|46|) is fulfilled. 


The instantaneous vacuum can be used as well for set¬ 
ting initial conditions. The instantaneous vacuum choice 
relies on the minimal energy state of the system defined 
by the Hamiltonian. Therefore, the requirement 

k ‘ rt = 0, (49) 

V =r hnit 

[55] for instance). In fact, one 
nlt such that both conditions 


(2fc etf )3/ 2 


has to be satisfied (see Ref. 
can find that there exists r/ 

( 46) and (491, are fulfilled. For reasons of comparability 


of the two approaches, we use this choice. In such a case, 
any difference between the two approaches is due to the 
higher order contribution to (48), which is present in case 


of the adiabatic vacuum-type normalization. 

The conditions for the validity of both the instantaneous 
and WKB vacua depend strongly on the evolution of the 
cosmological term z"/z during the pre-bounce contract¬ 
ing phase. A direct calculation leads to 

— = — a 2 ( m 2 — 2 H 2 + 2 nm 2 — + - /tf7 fi 2 
z V H 2 


-K 2 f7 2 



(50) 


This expression is valid at all times. In order to ana¬ 
lyze the shape of the effective potential it is convenient 
to divide the evolution into three background phases as 
mentioned in Sec. III. Then, analytical approximations 
for every phase can be used respectively. During the pre¬ 
bounce classical contracting phase, when p(t) p c , the 
scalar field undergoes an oscillatory behavior with an am¬ 
plitude proportional to \Jp(t). The Hubble parameter H 
is proportional to ^fp as well, whereas 17 ~ 1. Inserting 
these solutions into Eq. (50) yield terms which are pro¬ 


portional to different orders of ^fp. Averaging over the 
oscillatory contributions, which all have a characteristic 
oscillation time of 1/m, gives 


= —a 


(jn 2 — aiy/Hmy/ p(t) + a 2 Kp(i) 


a 3 np c 



(51) 

































where the constants a.i are determined by the averaged 
oscillations. Since the energy density is increasing for all 
times in the remote past, p(t) becomes sufficiently small 
in the remote past. Thus the m 2 -term will dominate for 
early times and therefore z"/z oc —m 2 a 2 . On a loga¬ 
rithmic scale as a function of ln(a/a B ) like in Fig. [2] the 
absolute value of the effective potential is then given by 
a straight line with gradient —2 and with a h\(\z"/z\)~ 
intercept of 21n(a B m) = —27.26. This result is obtained 
as well by a purely analytical analysis which is presented 
in Fig. [2] We use that the Hubble parameter H is ap¬ 
proximated by H(t) = for the pre-bounce 

phase when neglecting the fast oscillations, where H 0 de¬ 
notes the initial Hubble parameter. H 0 is determined by 
the mass and the parameter ck, namely H 0 = — m/3a. 
Integration leads to the analytical solution of the scale 
factor in the pre-bounce phase 

2 

a(t) = a* ^2 - — tj 3 , (52) 

where the prefactor a* is the scale factor for t = 
a/m. With this expression the analytical solution for 
Yn\(z"/z)\ = In a(t) 2 m 2 reads on the logarithmic scale 


In 


„, 4, / m 

= 2 In a* -I— In ( 2- 1 

3 V a 


2 ln?n. 


(53) 


As a function of ln(a/a B ), one gets the red line on the left 
in Fig. [2] This analytic solution is valid until the energy 
density starts to dominate over the constant mass term 
in Eq. ( |51| ) . The first term which is comparable to to 2 is 
proportional to v fp. With the analytic solution of yfp in 
the pre-bounce phase, 


\ pit) 


Ta 


-1 


1 — A-. [Yrrf + | sin(277if + 20 o )] 


(54) 


we can compute the time when the ^/p-term crosses ‘m 2 ’ 
in its amplitude. The oscillation term in Eq. (54) is av¬ 


eraged over T = 1/m. This transition point is referred 
to N pre := ln(a pre /a B ) = —5.38 in the figure. 

During the bouncing phase the potential energy param¬ 
eter x , is very small compared to the kinetic potential 
parameter j/, since we consider a kinetic bounce scenario. 
Then, in particular, x 2 <C y 2 is satisfied during this phase 
and the Hubble parameter can be reduced to 


H 2 = fy 2 (1 - y 2 ) ■ 


(55) 


The analytic solution for y around the bounce is given by 
y{t) = (1+3K/9 c (f—t B ) 2 ) -1 / 2 , as presented in Ref. [33] and 
the scale factor is related to y via a = a B \y\~ 1 ^. With 
these approximations the expression for z"/z reduces to 


z 


u 


z 


A- 1 

np c ai V 

(*) 

| 2 + 23 

(i 


(a) 

0 

3 ( 

jlV 

a b ) 

3 -1 


Note that this expression is positively valued and di¬ 
verges at the bounce. The absolute value of this ex¬ 
pression on the logarithmic scale and as a function of 
ln(a/a B ) provides the red lines around the bounce in 
Fig. i These approximations are valid until x 2 becomes 
significant in comparison to y 2 , let’s say x 2 > y 2 /10. 
The analytic solutions for x and y provide these tran¬ 
sition points of validity respectively before and after 
the bounce, namely N pTeB := ln(a preB /a B ) = —4.41 and 
N pmtB := In (a pmtB /a B ) = 3.62, as shown in Fig. [2] 

During slow-roll inflation, H ~ 1 such that z'^fz takes 
its classical expression. This leads to z" / z = (2 + 6e H — 
36 H )/r] 2 , with e H and <5 H the first and second Hubble flow 
functions, both much smaller than unity during infla¬ 
tion. Furthermore, a a l/») and therefore In [z"/z) oc 
21n(a/a B ), see Fig. [2j In particular the effective poten¬ 
tial is then given by z"/z = \a 2 H 2 . The curve can be 
approximated from the beginning of slow-roll inflation, 
i.e. when L = t B + f/m where / is an analytical ex¬ 
pression related to the Lambert function and t B is given 
analytically as well. For this time m = a B F _ 3 and the 
logarithm of the absolute value of the effective potential 
is given by \n\z"/z\ = ln(2/7y 2 ) = ln((l/2 )a 2 H 2 ). The 
value of the Hubble parameter H { is given analytically as 
well, see Ref. [33(. The approximation is valid starting 
from N post := ln(a post /a B ) = 5.53. For the approximation 
we use that 


H{t) = H, 


, e r . 

1- m(t — ti) 

x , 


(56) 


during slow-roll inflation, where e is the sign of the cosine 
of the phase parameter between the potential and kinetic 
energy parameters at the transition point between the 
pre-bounce and bouncing phases, and x t is the value of 
the potential energy parameter at t { . Furthermore, the 
scale factor undergoes an exponential growth with coor¬ 
dinate time, namely 


a(t) =a 1 e _ ^ (t - tl)(erm(t - t,) - 2xi) . 


(57) 


The analytic fit given by these two functions is displayed 
by the red line on the right side in Fig .[2 as a func¬ 
tion of ln(a/a B ) with a provided by Eq. (57). The slight 
difference between numerical results and the analytical 
solution during slow-roll inflation is due to the fact that 
the analytical approximations for this phase goes back on 
approximations even for the pre-bounce phase. Hence, 
small differences are propagated. 
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FIG. 2. Evolution of the cosmological term z" fz as a function 
of the number of e-folds ln(a/a B ), with m = 1.2 x 10~ 6 m P i. 
The parameters for the background are set as in Fig. [T] Dur¬ 
ing the pre-bounce contracting phase and slow-roll inflation 
In (z"/z) oc ±21n(a/a B ) + const. 


VI. THE SCALAR POWER SPECTRUM 


For scalar modes, the primordial power spectrum at 
the end of inflation is defined in terms of the mode func¬ 


tions by virtue of the definition (38). As we shall see 


in the next section, three ranges of wavenumbers can be 
identified, depending on how they compare to the effec¬ 
tive potential z" / z: (i) the infrared regime, (ii) interme¬ 
diate scales and (iii) the ultraviolet regime. 


A. The infrared regime 


The infrared limit (IR) of the primordial power spec¬ 
trum corresponds to modes such that k 2 <C \z" / z\ dur¬ 
ing the pre-bounce contracting phase. These modes are 
frozen during the bouncing phase and slow-roll inflation. 
The transition between the contracting phase and the 
bouncing phase occurs when H ~ —m/3, as discussed in 
Ref. (33j. At the transition, the effective potential term 

z" / z is well approximated by a"/a, since <j>/H remains 
nearly constant. This allows us to introduce the following 
IR scale (see Ref. j55|): 


kjTt :— 


Zy/2 


^m 2 


1/3 


4.7 x 10 _5 m 


(58) 


where the numerical value has been obtained for m = 
1.2 x 10~ 6 m P i, a B = 1 and p c = 0.41mp r 
There is an important difference with respect to the case 
of tensor modes. For small enough values of k in the 
IR regime, the tensor power spectrum tends to be scale- 
invariant. This is due to the fact that initial conditions 
for tensor modes are set when all the modes of inter¬ 
est are sub-Hubble (or, more precisely, k 2 /g> a"/a). For 
the scalar modes, however, it is impossible to set ap¬ 
propriate initial conditions at a time when all relevant 


modes are such that k 2 \z" / z\. Indeed, the condi¬ 
tions discussed in the previous section, (46) and (49), 
have to be satisfied as well respectively for the WKB and 
the instantaneous vacuum-type normalizations. In addi¬ 
tion, the absolute value of the effective potential term 
keeps decreasing in the remote past of the contracting 
branch. It is possible to find a time r/ Ilit , in the clas¬ 
sical contracting phase, at which (i) the absolute value 
of the effective potential term z" /z is close to zero and 
(ii) the conditions of validity of the vacuum states are 
fulfilled. Nevertheless for the WKB vacuum, at this par¬ 
ticular time, when \z" / z\ is minimal and condition (46) is 
satisfied, the effective potential does not strictly vanish, 
it is \z" / z\ = 2.1 x 10 _ 'Topi, and therefore only the modes 
with k > 4.5 x 10 _4 m P1 satisfy the condition k 2 > \z" / z\. 


B. The ultraviolet regime 


In the deformed algebra approach, the ultraviolet 
modes (UV) experience an exponential growth with in¬ 
creasing wavenumbers. This is due to the fl-factor in 
front of the wavenumber in Eq. ( |32| ) , which becomes neg¬ 
ative near the bounce. When approaching the bounce, 
the friction term in Eq. (33), namely H/H = l/(t — t B ), 
diverges. However, the approximate solution given in sec¬ 
tion IV shows that Uk vanishes at the bounce, since its 
generic expression is given by 

Uk = {t- t B ) [cich (k(t - t B )) - c 2 sh (k(t - t B ))], (59) 

where c\ and c 2 are numerical constants that have to 
be chosen in accordance to the initial state of the per¬ 
turbations. Thus, the equation of motion (33) has no 


singularity, and (H/H)iik oc y/k. We are left with the 
differential equation of an harmonic oscillator, but with 
an imaginary frequency and a constant term, say j3y/k, 
such that ilk + PVk — k 2 Uk = 0. Close to the bounce 
the generic solution to this equation is 


Uk — Pk ^ -f- cr_|_e^ -(- o_e 


— kt 


So, in the large k limit and close to the bounce the 
amplitude of scalar modes receives a real exponential 
contribution, which marks the Euclidean nature of the 
bounce (f2 < 0). For large scales however, one has 
|fife 2 1 -C \z" / z\ around the bounce. The solutions for 
the mode functions are 7 Ikiv) ~ + -B k f n drj'/z 2 (r]'), 

and the large-scale modes are thus qualitatively not af¬ 
fected by the Euclidean nature of the bounce. A similar 
behavior for the tensor modes was already discussed in 
Refs. HU and [333 - 

The characteristic energy scale /c uv beyond which the ef¬ 
fect of the Euclidean nature of the bounce qualitatively 
affects the evolution of the modes can be determined from 
an analysis of k 2 „ in Eq. ( f39| ). In vicinity of the bounce 
H — 1 and z"jz > 0. Therefore, a given mode has 
an imaginary time-dependent wavenumber for a certain 
period around the bounce, i.e. ‘k 2 {{ < O’. This is what 
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we call the ‘Euclidean phase’ in this approach. How¬ 
ever, the imaginary effective wavenumber only plays a 
role, if the interval of conformal time spend by the mode 
in this regime is large enough. Of course, in the Eu¬ 
clidean phase, it makes physically no sense to talk about 
time, the evolution parameter 77 however remains and 
quantifies the ‘period’ of the mode. If the mode spends 
more than one period in the region with complex effec¬ 
tive wavenumber, Ay rr < 0, the mode will be amplified 
significantly. This is in accordance with the analytical so¬ 
lutions to the approximated differential equation around 
the bounce, (331, which show an hyperbolic behavior for 
k(t — ts) 1. This leads to the following condition for 
the energy scale fc uv ; 


k uv Aij(k vv ) « 1. (60) 


A direct analytical analysis of this condition gives the 
following expression 


^uv — O’ 



V2 

2x72-1 


2.3 m P1 , (61) 


bounce. For such modes the effective description based 
on the continuous equations of motion might not be re¬ 
liable. Second, the observed amplification is due to an 
instability related to the elliptic type of the equation of 
motion for perturbations in the Euclidean regime. The 
Cauchy initial value problem might not be valid for the 
modes with k > fc uv which are strongly affected by the 
Euclidean nature of the deep quantum regime. 



k [ot p1 ] 


where the numerical value is obtained with use of a B = 1 
and p c = 0.41 ?n pl . 


C. Numerical Results 


FIG. 3. Primordial power spectrum for scalar modes in the 
deformed algebra approach for m = 1.2 x 10 _6 m P i and the adi¬ 
abatic vacuum for initial conditions in the pre-bounce (classi¬ 
cal) contracting phase. The cosmological background is fixed 
such that x B = 3.55 x 10“° and a B = 1 at the bounce. 


The scalar power spectrum is obtained by numerical 
integration of the equation of motion for the mode func¬ 
tions, respectively for the variable h and TZ for different 
phases in the time evolution, and the solution of the back¬ 
ground equations (12). Initial conditions for the pertur¬ 
bations can be set according to the WKB approximation 
referring to the adiabatic vacuum, or with the instanta¬ 
neous vacuum as shown in Sec. V. The initial conditions 
for the cosmological background are set in the contract¬ 
ing phase, such that the preferred value of the potential 
energy parameter x at the bounce is obtained. Note that 
the dynamics of the background and subsequently the 
shape of the power spectrum V s (k) are determined by 
the mass m of the scalar held, the value of the critical 
energy density p c and the phase 8q. 

Our numerical results are shown in Fig. [3] and Fig. [4] 
which display the primordial power spectra for the scalar 
modes, choosing the adiabatic (WKB) vacuum as initial 
conditions (Fig. [ 3 ]) and the instantaneous vacuum as ini¬ 
tial conditions (Fig. El). The three regions mentioned in 
the previous section (fc < fc IR , fc IR < k < fc uv , k > k uv ) 
can be well identified in the spectra. 

In the intermediate region (fc IR < k < fc uv ), the spectrum 
follows a characteristic oscillating behavior observed also 
in case of the tensor modes (see Ref. [33]). For the values 
of k > fc uv , the power spectrum is characterized by the 
exponential growth. This behavior should, however, be 
considered with care. First, the UV regime (k > fc uv ) cor¬ 
responds to the modes which are trans-Planckian at the 



FIG. 4. Primordial power spectrum for scalar modes in the 
deformed algebra approach for m = 1.2 x 10~ 6 m P i and the 
instantaneous vacuum for initial conditions in the pre-bounce 
(classical) contracting phase. The cosmological background 
is fixed such that i B = 3.55 x 10 -6 and a B = 1 at the bounce. 


The spectra for the adiabatic vacuum and the instan¬ 
taneous vacuum choices are almost identical. The only 
difference is a slight enhancement in the IR region for the 
instantaneous vacuum-type normalization, in comparison 
to the adiabatic vacuum choice. This effect is due to the 
difference in t’(.(binit) for both types of initial conditions. 
At values of fc < fc IR , the shape of the spectra is mostly 
due to the initial super-Hubble vacuum normalization. 
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In this limit ~ \ k 2 — — sa const, is almost the same 


for every mode at this particular initial time, because 
k 2 <C . Therefore, V s (k) = ^2 2 k^ t z 2 ^ ^ 3 - T ^ e 


initial V s (k) oc k 3 behavior is preserved in the further 
evolution due to the super-Hubble nature of the modes. 
Only the absolute amplitude changes, which is a result 
of the time dependence of the z parameter. 

Furthermore, it is worth stressing that at scales k < k uv 
the power spectra computed in this paper qualitatively 
agree with those obtained in the so-called “dressed met¬ 
ric” approach to perturbations in LQC [2S] ■ In that case 
the fl-factor does not appear in front of the Laplace oper¬ 
ator and the instabilities related to the Euclidean phase 
do not arise. Therefore, the corresponding spectrum at 
k > k vv becomes nearly-scale invariant in the “dressed 
metric” approach, as in the standard inflationary picture. 

We live the detailed associated phenomenology for a fu¬ 
ture study m- It is however important to stress that the 
spectrum derived in this article is basically in disagree¬ 
ment with data that are very precise for scalar modes. 
This rules out neither LQC in general nor the deformed 
algebra as a whole. But this shows that the set of spe¬ 
cific hypotheses presented here is obviously in tension 
with measurements. Ruling out a given setting is useful 
for future quantum gravity investigations. 


VII. DISCUSSION AND CONCLUSION 

In this work, the primordial scalar power spectrum in 
the so-called “deformed algebra” approach to perturba¬ 
tions in loop quantum cosmology has been derived. Our 
considerations were focused on the model with a massive 
scalar field. The instantaneous and adiabatic vacuum- 
type initial conditions were imposed in the contracting 
phase. The non-trivial issue in the evolution of modes is 
their behavior in the Euclidean phase (fi < 0) surround¬ 
ing the bounce. In this region, the equation of motion 
for the mode functions changes its type from hyperbolic 
to elliptic. In such a case, no preferred time direction 
exists. It is usually argued that a signature change im¬ 
plies instabilities because the oscillating time dependence 
‘exp(±iwf)’ is replaced by an exponential ‘exp(±wt)’ one. 
The growing mode leads to an instability if initial val¬ 
ues are chosen at some fixed t. The approach suggested 
in Ref. [13: precisely investigates how a boundary value 
problem for this kind of elliptic equations can eliminate 
the instability. This is, certainly, a path worth investi¬ 
gating. This article is devoted to the other hypothesis: 
considering seriously this real argument in the exponen¬ 
tial function and investigate its physical consequences. 
Even if the problem is not posed in the usual way for 
the partial differential equation in that case, reliable pre¬ 
dictions can still be obtained for sufficiently low values 
of k in the Fourier space representation. More pre¬ 
cisely, a characteristic scale fc uv discriminates between 
the modes which are amplified due to the imaginary effec¬ 


tive wavenumber k el , around the bounce. This is a mathe¬ 
matical consequence of the equation of motion in Fourier 
space. Of course, the physical interpretation as having or 
not enough “time” to oscillate and feel the quantum ge¬ 
ometrical structure near the bounce is not possible with¬ 
out time. But this is not something fundamentally new or 
surprising in quantum cosmology/gravity. For the modes 
satisfying k > k vv , the elliptic nature of the equations 
becomes important, leading to an abnormal amplifica¬ 
tion of the power spectrum. The effect is the same as 
the one observed earlier in case of the tensor perturba¬ 
tions m- In turn, for k < fc uv , that is for large wave¬ 
lengths, the modes are not subject for a sufficiently long 
period to the negative effective potential, k 2 it < 0. Thus, 
the corresponding modes are not affected by the hyper¬ 
bolic amplification, as discussed above. In the regime, 
k m < k < fc uv , a typical oscillatory behavior is ob¬ 
served. In the IR limit, the shape of the spectrum is de¬ 
termined by the initial vacuum normalization and scales 
as V s (k) oc k 3 . This behavior is very different from the 
one observed in case of the tensor modes (see Ref. 133]), 
where the power spectrum becomes nearly-scale invari¬ 
ant while k —>■ 0. This is because the massive scalar field, 
oscillating in the contracting branch, effectively behaves 
as dust matter. As it is known, the freezing of massless 
modes during such an evolution leads to scale-invariance 
of the power spectrum, as for the case of tensor pertur¬ 
bations. For the scalar perturbations in a model with 
a massive scalar field the gauge-invariant degree of free¬ 
dom v s is explicitly massive leading to a breakdown of 
the scale-invariance. 

Several points of the picture presented in this study still 
need to be addressed. First, the observational conse¬ 
quences of this calculations should be studied into the 
details. The key point for phenomenology is the knowl¬ 
edge of full duration of inflation. This is what translates 
coordinate wavenumbers used in this study into physi¬ 
cal wavenumbers that can be compared with data. The 
number of inflationnary e-folds is in one-to-one corre¬ 
spondence with the value of the scalar field at the bounce 
which, itself, depends on the phase of the oscillations in 
the contracting branch. The higher the field value (and 
therefore the fraction of potential energy) at the bounce, 
the longer the inflation period and the smaller the physi¬ 
cal scales at the bounce that are nowadays probed by as¬ 
tronomical observations. We leave a detailed phenomeno¬ 
logical study for a future work, 132 , but it is easy to 
guess from the primordial power spectrum derived in this 
study that the model as it is here investigated disagrees 
with data. It could simply be a consequence of the way 
modes with k > 1 are handled. But whatever the reason 
this problem should be stressed. Second, other propos¬ 
als for setting initial conditions should also be consid¬ 
ered. Here, the subtle issue of the very meaning of time 
in the Euclidean phase were deliberately ignored: modes 
were naively propagated through the Euclidean phase. A 
proper addressing of the well-posedness is crucial to ob¬ 
tain stable solutions in the k > fc uv regime (even if their 
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physical meaning is not clear due to the breakdown of 
validity of the effective equations under considerations) 
[35]. Furthermore, the matter content considered in this 
paper is no more favored by the observations of the cos¬ 
mic microwave background radiation. A caraful analy¬ 
sis of different inflationary potentials would therefore be 
desirable. In particular, the Colemen-Weinberg poten¬ 
tial with an unstable state may lead to inflationary spec¬ 
tra being in agreement with the up-to-date observational 
data. Such a change of the potential function would un¬ 


avoidably affect our predictions regarding the shape of 
the power spectra. 
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